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Abstract 



The deuteron is described as a quantum state on a ten-dimensional manifold M\q of Skyrme fields of 
degree two, which are obtained by calculating the holonomy of SU(2) instantons. The manifold M\q 
includes both toroidal configurations of minimal energy and configurations which are approximately 
the product of two Skyrmions in the most attractive relative orientation. The quantum Hamiltonian 
is of the form — A + V, where A is the covariant Laplace operator on M\q and V is the potential 
which Mio inherits from the Skyrme potential energy functional. Quantum states are complex- 
valued functions on the double cover of M\q satisfying certain constraints. There is a unique bound 
state with the quantum numbers of the deuteron, and its binding energy is approximately 6 MeV. 
Some of the deuteron's electrostatic and magnetostatic properties are also calculated and compared 
with experiment. 
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1 Introduction 

A fundamental challenge in particle physics is to understand nuclear forces from first princi- 
ples. So far, there is no quantitative understanding of nuclear binding starting with quarks 
and QCD, but one may attempt to investigate nuclei in terms of an effective low energy 
theory, like the Skyrme model (![]. In the Skyrme model, nucleons are solitons, and the 
parameters of the model are fixed so that the masses of the nucleon and the delta resonance 
are in agreement with experiment. All of nuclear physics can in principle then be derived 
from the Skyrme model (assuming that it is at least approximately right). 

Further simplifications are necessary to make progress with this programme. It is gen- 
erally agreed that treating the Skyrme model as a quantum field theory is very hard, and 
to study nucleons and their interactions one needs to reduce the Skyrme model to a finite- 
dimensional quantum mechanics. This can be done by picking out, as naturally as possible, 
a 6iV-dimensional set of Skyrme fields to model N nucleons. One nucleon is a quantised 
state of a Skyrmion, the lowest energy classical configuration with unit baryon number. The 
Skyrmion has six classical degrees of freedom - three translational and three rotational - and 
their quantisation gives a nucleon with momentum and correlated spin and isospin. There 
are various suggestions for a twelve- dimensional set M\2 of two-Skyrmion fields. All include 
the product of two well-separated single Skyrmions (as suggested by Skyrme in |fj]) to model 
well-separated nucleons, but it has been known for years that the product ansatz is not good 
for Skyrmions close together. The minimal energy two-Skyrmion fields are configurations 
with toroidal symmetry in which two separated Skyrmions have coalesced. These should be 
included in M 12 . 

The idea we favour, in principle, is to define M 12 as the unstable manifold of the hedgehog 
two-Skyrmion configuration. This is discussed in detail in |J. There is a quantum Hamilto- 
nian for the motion on M12, of the form H = —A + V, where A is the Laplacian on M12 and 
V is the potential energy. In this paper we do not consider this Hamiltonian, but a further 
simplification. It is known that Myi has two naturally defined, low-lying submanifolds, which 
should be most important for low-energy physics of two nucleons. The first, denoted Mg, is 
eight-dimensional and consists of the toroidal configurations of minimal energy. Braaten and 
Carson considered the quantum mechanics on M 8 in [Q. They calculated the Hamiltonian 
and found the lowest energy stationary state compatible with the nucleons being treated as 
fermions. This state is identified with the deuteron. Not surprisingly, Braaten and Carson 
found a very large binding energy, and a small physical size. This is because the Skyrmions 
are not allowed to explore the full twelve- dimensional M 12 , but only the eight-dimensional 
Mg. Nevertheless, Braaten and Carson's analysis is very instructive and is an important 
inspiration for our work. 
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We have considered the quantisation on a ten-dimensional submanifold Mio of M12. 
Again the Hamiltonian is derived from the kinetic and potential parts of the Skyrme model 
Lagrangian, restricted to the fields in M\q. The manifold Mio is the set of attractive channel 
fields, which is a low-lying valley in M 12 . It can be defined by taking two well-separated 
Skyrmions, oriented so that they maximally attract (i.e. one is rotated by 7r relative to the 
other, with the axis of rotation orthogonal to the line joining them), and allowing them to 
approach until they coalesce into a toroidal configuration. Ideally one would follow a gradient 
flow curve, or a path of steepest descent. The ten coordinates of Mio are accounted for by 
the separation parameter, overall translations, rotations and iso-rotations. Mio should be 
thought of as the smallest subset of M i2 which includes Skyrmion separation. 

In practice, we have not followed the prescription above for constructing M 10 . Instead, 
we have defined an approximation to M 10 consisting of Skyrme fields generated by Yang- 
Mills instantons. The details of the construction are explained in section 3. The instanton 
approach has some advantages which compensate for the fact that the fields do not have 
quite as low energies as those obtained by gradient flow. The main advantage is that the 
fields and their currents can be calculated by numerical integration of ordinary differential 
equations, whereas the gradient flow approach requires the numerical solution of a partial 
differential equation. Also, there is an explicit separation parameter p in the instanton data, 
which is a convenience. A disadvantage of the instanton approach is that it requires the 
pions to be regarded as massless. Other approaches can deal with pions which have their 
physical mass. It would certainly be worthwhile to work with M w as defined by the gradient 
flow, and compare with the results here. 

It is worth comparing the picture of the two nucleon interaction that emerges from the 
Skyrme model with conventional nuclear potential models. The basic difference is that 
no subset of the coordinates on M 10 can be identified simply with the positions of the 
two nucleons. While the two Skyrmions are separated one can identify two points where 
the energy density (or baryon number density) is maximal, but when they coalesce these 
points disappear as the Skyrmions lose their identities. The energy density of a toroidal 
configuration is maximal on a circle. Related to this is the fact that, in the attractive 
channel, two Skyrmions are never closer together than in a toroidal configuration, and we 
shall define below a separation parameter p lying in the range [po, oo], where po corresponds 
to the toroidal configurations, and is positive. When Skyrmions approach head-on they pass 
through a toroidal configuration, scatter through 90° and separate again. The radial part 
of the deuteron wavefunction u(p) satisfies a radial Schrodinger equation on the interval 
[pojOo). Since M w is smooth at po, the boundary condition there is that u is finite and 
du/dp = 0. In fact, u(p) is maximal at p , which is where not only the potential V, but 
even the effective potential V cS occurring in the radial equation is deepest. This is quite 
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different from a point-particle description of nucleons where a hard core potential is required 
to keep the nucleons from being close together. With a hard core, the boundary condition 
is it(0) = 0, or if the core potential is infinite for p < pi, then u(pi) = 0. In this situation, 
the wavefunction is usually maximal outside the point where the potential is deepest. As a 
consequence, a stronger attraction is needed to compensate. The Skyrme model seems to 
solve the usual problem of a 'medium range central attraction' through the geometry of Miq. 

As an aside, we point out that the quantisation of a ten-dimensional family of two- 
Skyrmions was also considered in ||. The ten-dimensional family, like our manifold Miq, 
was obtained by acting with translations, rotations and iso-rotations on a one-parameter 
family of fields which interpolates between the toroidal configuration and two-Skyrmions 
which are approximately of the product form, but have additional reflection symmetries. 
These reflection symmetries also play a role in our analysis, but our quantisation scheme 
differs from that of ||, which appears incorrect. In H an effective radial potential between 
Skyrmions was calculated on the interval [0, oo). The wavefunction was required to vanish 
at p = 0, and there was (by a long way) no bound state. We find a bound state with a very 
similar potential because of our different (and we believe correct) boundary conditions. 

There is well-defined angular momentum in our formalism, but because the Skyrmions 
are extended objects and interact, there is no well-defined split into orbital and spin parts 
of the angular momentum. Consequently we cannot perform the usual decomposition of the 
deuteron wavefunction into s-wave and d-wave parts. Nevertheless we can calculate physical 
quantities related to this splitting, like the electric quadrupole moment. 

The paper is organised as follows. Section 2 contains a review of the Skyrme model and its 
symmetries. In section 3 we review how Skyrme fields are obtained from instantons, and in 
particular how 2-Skyrmion attractive channel fields are obtained and parametrised. Section 
4 describes the topological structure and the symmetries of the manifold M 10 . The form of 
the Skyrme Lagrangian restricted to attractive channel fields is described in section 5. The 
Lagrangian depends on nine functions of the Skyrmion separation parameter p. These are 
calculated numerically, but for large values of p we also have analytic formulae, which provide 
a useful check. The appendix contains details of the numerical methods used. In section 
6, the quantum Hamiltonian in the attractive channel is derived from the Lagrangian, and 
it is explained, using angular momentum analysis, how the stationary Schrodinger equation 
reduces to a radial equation. The bound states of this radial equation are discussed in 
section 7, and the lowest physically allowed state is identified with the deuteron. In section 
8, various physical properties of our deuteron state are calculated and compared with the 
values obtained by Braaten and Carson, and with experiment. Some of the results are 
satisfactory, others less so. For the first time a sensible binding energy of order 5 MeV is 
obtained in a Skyrme model calculation. The deuteron has the right size, too, but its electric 



4 



quadrupole moment is too large and its magnetic dipole moment is too small. We conclude, 
in section 9, with some remarks on how one might correct our calculations on M\q without 
going to a full quantisation on the manifold Myi- Ideally, however, quantum mechanics on 
M 12 is what should be done next. 



2 The Skyrme Model 



We set the speed of light to 1 and use the metric diag( 1,-1,-, 1,-1) on Minkowski space. 
Points in Minkowski space are written as (t, x) with coordinates x M , fi = 0, 1,2,3, and the 
Einstein summation convention is used throughout. The basic field of the Skyrme model is 
the 577(2)-valued field U(t,x), which can also be expressed in terms of the pion fields: 

U(t,x) = (r(t,x)+iir(t,x)-T (2.1) 

where r l5 r 2 , r 3 are the Pauli matrices and a 2 + 7r 2 = 1. The equations of motion for U are 
the Euler-Lagrange equations derived from the Lagrangian density 

C = -^tr(R^) + -L^ tT ({R^R u }{R^R»}). (2.2) 

The right- currents R^ are defined via 

R„ = (d^U)U\ (2.3) 

where denotes partial differentiation with respect to x M . The constants F n and e are 
free parameters of the Skyrme model. We found it convenient to remove them from our 
calculations by using F^/Ae as the unit of energy and 2/eF^ as the unit of length. Thus we 
work in geometrical units, and the Lagrangian density takes the form 

C = -itr(i^) + ±tr([R p ,R 1/ }[R»,R>'}). (2.4) 

The parameters F n and e can be fixed in a number of ways. Most authors adopt the 
approach of |j) and |7| where F n and e are tuned to reproduce the masses of the proton and 
the delta resonance without and with the physical pion mass respectively. All papers written 
so far on the deuteron in the Skyrme model take into account the physical pion mass and 
use the set of parameters given in J7[. For ease of comparison we will follow that practice, 
although the pion mass is assumed to be zero here. Thus, with the values of F w and e as in 
pj, our units are related to conventional units via 

F 2 

-^ = 5.58MeV and — = 0.755 fm. (2.5) 
4e eF % K ' 



Since 



n = 197.3McY tin - Ki.s ( ^ ^-|_ ) . (2.6) 
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it follows that, in our units where F n /4e = 2/eF n = 1, h = 46.8. 

It is useful to think of the Skyrme model as an infinite-dimensional Lagrangian system 
whose configuration space Q is the space of maps 

U : R 3 i-> 577(2) (2.7) 



which obey 



hm U(x) = l 2 . (2.8) 

\Ju\ — ^OO 



The Lagrangian L = J d 3 xC has the usual form L = T — V , where the T is the kinetic energy 

T = J d 3 x{-^tr(iV?o) - ^tr({R h R }{Ri, R ])} (2.9) 
and V is the potential energy 

V = Jd 3 x{-~tr(R i Rj-^tr([R i ,R j ][R h R j })\. (2.10) 



The condition ( |2.8[ ) is imposed in order to ensure that the elements of Q have finite 
energy. It is important topologically because it allows one to compactify R 3 to S 3 and to 
regard an element U of Q as a map S" 3 i— > SU{2) = S 3 . Thus ir (Q) = n 3 (S 3 ) = Z, showing 
that Q is not connected but has components Qb labelled by an integer B which is the degree 
of any element of Qb. Physically, the integer B is interpreted as the baryon number of a 
Skyrme field. It can be calculated by integrating the zeroth component of the (conserved) 
baryon number current 

B^ = ——tx{R v R a Rp), (2.11) 
where we use the convention 60123 = — e 0123 = 1. Thus, 

B° = -B = --^triR^Rk) = - JLtr([i2!, R 2 }Rs) (2.12) 



B* = Bi = %£tT{RjR k Ro). (2.13) 



and 



Then the baryon number is 

B = J d 3 xB°. (2.14) 

The fundamental group of Q is also important here. In general, the fundamental group 
of a space depends on the choice of a base point, but it is explained in § that 

7Ti(Q) = 7Ti(Q B ) = Z 2 (2.15) 
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for any B G Z. The group 7Ti(Q) is generated by the rotation of a single Skyrmion by 2tt. 
More generally, it is shown in || that the loop generated by rotating an element oiQs by 2n 
is contractible if B is even and not contractible if B is odd. A further important topological 
result, due to Finkelstein and Rubinstein |§, is that the rotation of one Skyrmion by 2ir and 
the exchange of two Skyrmions are homotopic paths in Qb, for \B\ > 2. 

We will use the following conventions when describing Skyrme fields. We refer to the 
domain R 3 of a static Skyrme field as physical space and to the range SU (2) as iso-space. 
We write {e±, e%,e%\ for the canonical orthonormal basis of R 3 (ei = (1,0,0) etc.) and in 
the decomposition (|2.1|) of elements of 577(2) we refer to ir a (a = 1, 2, 3) as the a-th axis in 
iso-space. G stands for an SU(2) matrix and D(G) for the SO (3) matrix associated to G 
via 

D(G) ab = \tx{r a G n G^. (2.16) 

For fixed B, the symmetry group of the Lagrangian system with configuration space Qb 
and Lagrangian L is 

P x SO(3Y x E 3 . (2.17) 

Here P is the combined parity operation in space and iso-space, 

P : U(x) i-> U\-x), (2.18) 

and SO(3Y is the group of iso- rotations. Its action can be written in terms of an SU{2) 
matrix C as 

U{x) i — ► CU(x)Cfl (2.19) 

or, in terms of the pion fields, 

ir(x) i-> D(C)tt(x). (2.20) 

Finally, the euclidean group E 3 is the semidirect product of the spatial rotation group SO (3) J 
and the group R 3 of translations. An element (D(G), S) G E 3 acts on a vector x according 
to 

x^D(G)x + S (2.21) 

and on Skyrme fields via pull-back: 

U(x) h-> U{D(G)~ l (x - S)). (2.22) 



The group ( |2.17 ) has discrete subgroups which will be important for us. We therefore 



introduce the notation O a ; (a, % G 1,2,3) for simultaneous rotations by ir around the a-th 
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axis in iso-space and the i-th axis in physical space. It will also be useful to write Oq{ for 
the rotation by it about the i-th axis in physical space only and O a o for the rotation by 7r 
about the a-th axis in iso-space only. Finally we define hyperplane reflections in space and 
iso-space P ai = PO ai . 

For the discussion of the kinetic energy we need carefully to define angular velocities: the 
left- invariant or body-fixed angular velocity for spatial rotations oj is defined via 

w-t = G*G, (2.23) 

where U = —\Ti (so that [ti,tj] = e^tk) and the dot denotes differentiation with respect to 
time. Equivalently, in terms of D(G): 

u-s = D(G)- l D(G), (2.24) 

where the Sj are 3x3 matrices with components (sj)jfc = —£ijk- They also satisfy [si, Sj] = 
EijkSk- The corresponding right-invariant or space-fixed angular velocity is less important 
for us; it is given by D(G)uj. 

Similarly one defines the body-fixed angular velocity f2 for iso-rotations via 

ft-t = C^C. (2.25) 

The corresponding space-fixed angular velocity is again given by D(C)Q. The angular 
momenta associated to the angular velocities in space and iso-space will be discussed in 
section 6. 



3 Attractive Channel Skyrme Fields from Instantons 



The idea that Skyrme fields can be generated from Yang-Mills instantons was proposed in [10 



and subsequently developed in a number of papers, including JTI| and [fT2| . The attractive 



channel fields generated from instantons were also studied in . 

Briefly, consider the fourth component of a self-dual SU(2) Yang-Mills field, or in- 
stanton, on R 4 , and define the holonomy along all lines parallel to the X4-axis 

/oo 
A±(x, x 4 )dx4 (3.1) 
-oo 

(the right hand side is a path-ordered exponential). The formula (|3.1| ) defines an SU(2)- 
valued field which can be regarded as a Skyrme field on R 3 . Provided A± decays sufficiently 
rapidly as |cc| — ► oo, U(x) satisfies the asymptotic condition Q2.8|). Moreover, the degree, or 
baryon number of U equals the instanton number or charge [O] . 
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The expression ( |3.1| ) is formal. It is computed by solving the differential equation 

d 4 U = -A 4 U, (3.2) 
where U is defined on R 4 . One imposes the boundary condition 

lim U(x, X4) = 1 2 (3.3) 

X4 — > — OO 

and defines the Skyrme field U(x) via 

U(x) = lim Uix.xA. (3.4) 

For details of the numerical implementation of these steps we refer the reader to the appendix, 
and also [^3[ . 



The simplest Skyrme field one obtains this way (apart from the vacuum) is the hedgehog 
Skyrme field of degree 1. This is generated by a unit charge instanton centred at the origin. 



As shown in [[10) the equation (|3.2|) can be solved analytically in this case, and gives the 
Skyrme field 

U H (x) = expif(r)x-T (3.5) 

where 

/(r)=7rfl-(l + A)-^. (3.6) 

The scale parameter A is chosen to minimise the energy of the field, which is interpreted 
as the Skyrmion mass. The minimum is obtained for A = 2.11; Then the mass is M = 
1.243 ■ 127r 2 = 147.2 and the moment of inertia is A = 8.363 -167r/3 = 140.1 (normalised 
so that the kinetic energy of a hedgehog spinning with frequency u is \ku) 2 ). Using the 
formula M + 3fo 2 /8A ||, it follows that the prediction for the nucleon mass is M = 153.1, 
which is 854.1 MeV in physical units. This rather small value results from the combination 
of the instanton generated profile function ( |3.6| ) with the choice of units ( |2.5| ) designed to 
reproduce the experimental values for the nucleon and delta masses when the physical pion 
mass is taken into account. However, the absolute value of the nucleon mass (and later the 
deuteron mass) is not important here. We are mainly interested in the deuteron's binding 
energy, which is the difference between its own mass and the mass of two free nucleons. For 
later use we also note that / has the simple form 



m - £ (3.7) 



when r is large. The constant p (which may be interpreted as the Skyrmion's pion dipole 
strength) is 7rA/2. 



9 



The attractive channel two-Skyrmion fields, generated by instantons, were discussed in 



detail in ||12|| . There is also an analytic formula for A4, given by Jackiw, Nohl and Rebbi 
(JNR) in |T4[], but the equation (|3.2|), which gives the Skyrme fields, must still be solved 
numerically. The JNR formula is somewhat gauge dependent and ungeometrical, so before 
giving it we recall Hartshorne's geometrical description of charge two instantons. Hartshorne 
showed that associated to any charge two SU(2) instanton in R 4 there is a circle and an 
ellipse in R 4 (the circle can degenerate to a straight line but this will not be relevant here). 
The ellipse is in the same plane as the circle, and interior to it. Moreover, the ellipse satisfies 
the Poncelet condition; that is, there exists a triangle with vertices on the circle whose sides 
are tangent to the ellipse. Poncelet's theorem states that if one such triangle exists, then 
there are infinitely many and any point on the circle can be taken as one vertex. 

Two well-separated Skyrmions maximally attract if one is rotated relative to the other 
by 7r about a line perpendicular to the line joining them. The attractive channel interpolates 
between such Skyrmion pairs and the toroidal two-Skyrmion, where the Skyrmions coalesce 
into the minimal energy configuration. The attractive channel is modelled by a subset 
of the Skyrme fields generated by charge two instantons. This is the subset where the 
Hartshorne circle lies in a spatial plane in R 4 (orthogonal to the X4-axis) and where the 
ellipse is concentric with the circle. Without loss of generality, we may suppose the common 
centre is at the origin in R 4 , that the circle is in the Xiavplane, and that the ellipse has its 
major axis along the 1-axis, see figure 1. In this standard form there are two parameters in 
the Hartshorne data for attractive channel fields. These may be taken to be the radius L of 
the circle and the length L\ of the semi-major axis of the ellipse. We denote the length of 
the semi-minor axis by L 2 . The Poncelet condition is that L\ + L 2 = L. The configuration in 
figure 1 is rather symmetric, and this manifests itself in the symmetries of the corresponding 
Skyrme fields, which we will discuss in section 4. 

When Li 3> L 2 , the instanton associated with the circle and ellipse of figure 1 generates a 
Skyrme field consisting of two well-separated Skyrmions, centred near the ends of the ellipse 
(the ^-coordinate is irrelevant now). When Li = L 2 = |L, the ellipse becomes a circle, so 
there is higher symmetry, and the instanton generates a toroidal Skyrme field. 

To progress, we need the explicit JNR formula for A± for the instantons whose Hartshorne 
data is as in figure 1. To this end we choose one triangle of the Poncelet porism; following 
Hosaka et al. [[HJ we choose the isosceles triangle shown. The JNR data consists of the 
points Xi, X 2 , X 3 in R 4 as shown and associated weights Aj, A 2 , A3. These can be rescaled so 
that A2 = A3 = 1. It is natural to define the angle 1? = arcsin(l/(l + Ai)), whose geometrical 
significance is shown in figure 1. In JNR terms, the essential parameters of the instantons 
are the radius L of the circle and the angle 1?, which lies in the range (0, |]. The semi-major 
axis of the ellipse is then L\ — (1 — sin$)L. The JNR gauge potential in R 4 has a fourth 
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component of the form 

At = idilnu(x)Ti, (3.8) 
where v is the following function of ) e R . 



n=l 



\x - X n \ 



There is a slight technical complication in that the JNR gauge potentials do not decay 
sufficiently rapidly at large distances to extend smoothly from R 4 to S 4 ; it turns out that 
an extra factor of —1 must be attached to U in (3A), see also [Q. The Skyrme fields in the 
attractive channel are then given by (|3.2|) - (|3.4|). 



The Hartshorne data in the standard form described above have two parameters, L and 
and acting with spatial translations, rotations and iso-rotations on the Skyrme fields 
generated from the corresponding instantons one obtains an eleven- dimensional family of 
Skyrme fields, which is not desired. Instead of retaining L and d as independent parameters 
we therefore fix L, for each value of to the value which minimises the Skyrme potential 
energy. This fixes the scale of the Skyrmions in the configuration, so that for *& « we 
obtain the product of two well-separated hedgehog Skyrmions with the scale factor A = 2.11 
as in ( |3.6|) and oriented so that the attraction is maximal. In the limit i? — > 0, the energy 
of this configuration tends to 2M = 294.4, twice the energy of a hedgehog Skyrmion. For 
•d = 7r/6 we obtain a toroidal Skyrme field whose energy is 1.19-247T 2 = 282.0. 

As a separation parameter we shall use p = 2L(1 — sin$) = 2L X . This makes physical 
sense: p is positive when i? = tt/6, taking the value po = 1.71, and it tends to infinity as 
i? — > 0. For small p is an accurate estimate of the distance between the Skyrmion centres 
(for example, 2L is not so good). The field approaches a product of hedgehog fields with this 
separation, as can be seen by considering another triangle in the Poncelet porism, namely 
one which is isosceles and has one vertex on the 1-axis (see also ||12|| ). Mathematically, our 
definition of separation is a little awkward, as L is a function of $ which is determined 
numerically, but this causes few problems in practice. 

To sum up, starting with the standard Hartshorne data depicted in figure 1 we generate 
a 1-parameter family of attractive channel Skyrme fields in standard orientation, which 
we denote by U(p,x). The iso-orientation of U(p,x) is chosen so that for large p it is 
approximately a product of two well-separated hedgehog Skyrmions on the 1-axis with one 
hedgehog rotated relative to the other by 7r about the 3-axis: 

U(p,x) w U H (x + ipei)r 3 U H {x - |pei)r 3 . (3.10) 

When p = p , U(p,x) is the toroidal Skyrme field, centred at the origin and with the axis 
of the torus along the 3-axis. 
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4 Topology and Symmetries of the Manifold of Attractive 
Channel Fields 



We define the manifold Mio of attractive channel fields as the orbit of the family of standard 
fields U(p,x) under the symmetry group (|2.17| ). For most of this paper we will not be 



concerned with overall translations of the standard field U(p,x). Thus we also define the 
manifold Mj of centred attractive channel fields as the quotient of M\q by the translation 
group R 3 . Alternatively one can think of M7 as the orbit of the family of standard fields 
U (p, x) under 

g = 50(3) J x SO(3) J x P. (4.1) 

For fixed p, the orbit is diffeomorphic to the quotient of Q by the isotropy group of its action 
on U (p, x). 

The isotropy group consists, by definition, of all elements of Q which leave U(p, x) invari- 
ant. The standard Hartshorne data in figure 1 is clearly invariant under reflections in the 
X1X2, x 2 x 3 and X1X3 planes, but the Skyrme fields U(p,x) are only invariant under certain 
combinations of these reflections with reflections in iso-space. Explicitly these combinations 
are: 

P21 ■ (tti, tt 2 , tt 3 ) h-> (tti, -tt 2 , tt 3 ) and (x 1 ,x 2 , x 3 ) i-> (-x 1 ,x 2 , x 3 ) 
P 22 : (vri, vr 2 , tt 3 ) (tti, -tt 2 , tt 3 ) and (x 1 ,x 2 , x 3 ) 1 ^ (xi, -x 2 , x 3 ) 
P33 ■ (tti, tt 2 , tt 3 ) i-> (7ri, tt 2 , -7T 3 ) and (xi, x 2 , x 3 ) ^ (an, x 2 , -x 3 ). (4.2) 

The group generated by these maps is an abelian subgroup of ( |2.17| ) of order 8 which we 
denote by Q$. Its elements are 

Gs = {1, O n , 12 , O 03 , PO 30 , P 21 , P 22 , P 33 }. (4.3) 

It is easy to check that 

Vier = {l, O n , O 12 ,O 03 } (4.4) 

is a subgroup which is isomorphic to the viergruppe. Using 

P0 3 qP 2 \ = On 
P0 3 qP 22 = 0\ 2 

PO 30 P 33 = O 03 (4.5) 
one shows further that there is an isomorphism 

G 8 Vier x Z 2 , (4.6) 
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where Z2 ={1, —1}, which identifies -PO30 G Gg with (1,-1) e Vier x Z2. It follows that for 
p > Po the isotropy group of C/ (p, x) is C/g and that the orbit under Q is diffeomorphic to 

(SO(3Y x SO(3) J ) /Vier. (4.7) 

When p = po, there is an additional invariance. The field U(p ,x) is invariant under 
spatial rotations about the 3-axis by some angle x an d simultaneous iso-rotation about the 
3-axis by 2%. Explicitly 

e- tXT3 U(p , D(e^ T3 )x)e lXT3 = U(p , x), x G [0, 2tt). (4.8) 

Such transformations form an 5*0(2) subgroup of SOfi) 1 x S'0(3) J which we denote by Ti. 
It follows that the isotropy group of U(po, x) is the semi-direct product of H with the group 
consisting of {1, P 2 2, P33, On}. Thus the orbit under Q is diffeomorphic to 

{SOiZ) 1 x SO(3) J ) /0{2). (4.9) 

Here 0(2) is the semi-direct product of 7i with the group {1, On}. Note that the orbit under 
the action of Q is six-dimensional if p > p but only five- dimensional if p = p . 

Concretely, any field in M w can be written in terms of SU(2) matrices C and G and a 
translation vector S as 

CU{p,D(G)-\x- 5))C f . (4.10) 

For p > po, D(C), D(G) and 5 can be used to parametrise the fields, but we should identify 
points related by the right action of O 3 and On: 

(D(C),D(G)) ~ (D(C),D(G)O m ) 

(D(C),D(G)) ~ (D(C)O 10 ,D(G)O i). (4.11) 

It follows, as a consequence of the group structure, that points related by the right action of 
O12 are also identified. For p = p we identify points related by the right action On and by 
the right action of elements in 7i: 

(D(C),D(G)) ~ (D(C)O w ,D(G)O 01 ) 

(D(C),D(G)) ~ ( J D(Oe- ixr3 ), J D(Oe-^ T3 )). (4.12) 

The ^-orbit structure of M 7 can be used to coordinatise M 7 in terms of the separation 
parameter p and Euler angles ((f), 6, 1^) for 50(3) J and (a, /5, 7) for .SO^) 7 . More precisely 
we write an element D(G) of SO (3) J as 

£1(0 = e ^ 3 e 0S2 e^ 3 (4.13) 
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where s« are the generators of £0(3) defined in section 2. Similarly an element D(C) of 
SO (3) 1 can be written as 

D(C) = e aS3 e^ 2 e^ 3 . (4.14) 

The range of the Euler angles is 0, i/;, a, 7 G [0, 27r) and 9, (3 G [0, 7r). However, points related 
by the maps O03, On, O12 should be identified. These maps now take the form 



O03: 




-A 




ana, 


7 1— > 


7, 






9 1- 


+ 0, 






1p 1— ► 


^ + 


7T 


On : 


/?•- 


-> 7T — 




a 1— > a + 7r, 


7 1— > 


-7, 








-> 7T — 




h-> + 7T, 


l/) 1 — ► 






O12 : 




-> 7T — 




a i— ► a + 7r, 


7 1— > 


-7, 






H 


-> 7T — 


0, 


h-> (f) + 7T, 




7T — 





(4.15) 



Using the orbit structure of M7 it is not difficult to understand its homotopy structure. 
This will be important in the quantum theory. All the results that are required here were 
already derived in ||, but it is useful to briefly state and prove them using our notation. 
Since M 7 is connected it is sufficient to consider paths which begin and end at the toroidal 
configuration U(p ,x). 

Our first claim is that loops generated by either spatial or iso-spatial rotations of U(p, x) 
by 27r are contractible in M7. To see this we first note that rotations by 2ir about any two 
axes are homotopic to each other; thus we can specify a convenient axis without loss of 
generality. The claim is clearly true for a loop generated by a spatial rotation of U(po, x) by 
27r about the 3-axis, since such a rotation is equivalent to an iso-rotation by Att about the 
3-axis. Since rotations by 47r are homotopic to the identity the loop is contractible. To show 
the contractibility of the loop generated by an iso-rotation by 2tt we also initially consider a 
rotation about the 3-axis. Then we smoothly increase p until U(p, x) is approximately the 
product of two well-separated hedgehog fields as in ( |3.1(J| ). Since the hedgehogs are rotated 
relative to each other about the 3-axis, iso-rotating the whole configuration by 2tt about 
the 3-axis is the same as iso-rotating each of the hedgehogs by the same amount about the 
3-axis. However, by our earlier remark this loop is homotopic to rotating each of hedgehogs 
by 27r about the 1-axis. Moreover, for hedgehogs iso-spatial rotations are the same as spatial 
rotations. Then, since the two hedgehogs are separated along the 1-axis, a spatial rotation 
of each by 2tt about the 1-axis is equivalent to an overall spatial rotation by 2tt about the 
1-axis. The contractibility of the loop then follows from our earlier result that a spatial 
rotation by 2ti is contractible. 

Our second claim is that all loops in M7 are either contractible or homotopic to a rotation 
of one of the hedgehog fields in the product ansatz ( ft.lO| ). It then follows that 7r 1 (M 7 ) = Z 2 , 
so M 7 has the same fundamental group as the space Q2 of all Skyrme fields of degree 2. To 
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prove the claim we have to consider the paths which connect points identified by O03, On 
and 0\2- First consider the path To3 in M7 obtained by performing a spatial rotation of 
U(po, x) by 7r about the 3-axis. The map O03 identifies the endpoint with the starting point, 
so the path is closed. However, for toroidal configurations a spatial rotation by it about the 
3-axis is equivalent to an iso-rotation by 2% about the 3-axis. Hence the path is contractible 
by our first claim. The second path we consider is obtained by rotating U(po, x) in both 
space and iso-space by it about the 1-axes. This path, denoted Tn, is closed because the 
map On identifies the endpoint with the starting point. By increasing p so that U(p, x) is 
of the product form ( [3.1 OP this loop can be deformed into the following 

e~^ x ' n U H {D{e^)x + \pe l ))nU H {D{e^)x - \pe x )r^ T1 
= U H (x + \pe 1 ))T 3 U H (D(e i ^)x - ^pe^, X e [0, tt). (4.16) 

Hence it is homotopic to a rotation by 2rr of one of the two Skyrmions in the product ansatz. 
We know that such a loop is not contractible in Q2 and, a fortiori, it is not contractible in 
Mf. It follows also that the loop T\ x is homotopic to a rotation by 47r of one of the Skyrmions 
and hence contractible. Finally there is the loop T12 associated with the identification 0\i- 
It is obtained by performing a spatial rotation by i\ about the 2-axis and an iso-rotation by 
7r about the 1-axis. It follows from 12 = 0nOo3 that = TnTf^. Hence Ti 2 is in the 
same homotopy class as Tn and not contractible. 

To summarise, loops in M7 are either contractible or homotopic to a rotation of one 
Skyrmion by 2rr. Our discussion of the loops Tn, r 3 and Ti 2 suggests the following physical 
interpretation of the identifications On,Oo3 and O12 in the asymptotic region of M7 which 
describes well-separated Skyrmions: On identifies configurations related by rotating one 
of the Skyrmions by 2tt, O12 identifies configurations related by the exchange of the two 
Skyrmions and O03 identifies configurations related by the exchange of the two Skyrmions 
and the simultaneous rotation of each by 7r. As mentioned in section 2, the rotation of one 
Skyrmion and the exchange of two Skyrmions are homotopic paths in Q 2 . We have seen 
that the two operations are also homotopic in M7. Note, however, that the homotopy does 
not require the creation of a Skyrmion - anti-Skyrmion pair, as is often claimed. In fact, we 
have checked that the baryon density for all fields in M7 is non-negative everywhere. 

5 The Lagrangian in the Attractive Channel 

The Lagrangian in the attractive channel is, by definition, the restriction of the Skyrme 
Lagrangian L to M\q. We call it L att . The evaluation of potential energy V gives a function 
of p only which we also denote by V. To calculate the kinetic energy, denoted T att , we allow 
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C, G, S and p in (|4.10 ) to vary with time. We write 

Ri = {diU)U\ R p = (d p U)U j and y = D(G)-\x - 5), (5.1) 

where d p denotes differentiation with respect to p; we also introduce the notation R for the 
vector with components R^. We then find, in terms of the angular velocities u> and Q defined 
in section 2, 

R = C ([-fft-r, U]fr -u-yxRy- S-D(G)R y + pR p ) C ] . (5.2) 



Here, all functions are evaluated at y and we have written R y to indicate that the differen- 
tiation should be carried out with respect to y. Inserting this formula into ( |2.9| ), changing 
integration variables to y (d 3 x = d 3 y) and then changing the name of the integration variable 
from y back to x, we find that the kinetic energy can be calculated solely from the standard 
field U(p, x) and its currents Ri and R p . 

The calculation is simplified by symmetries. It follows from the invariance of the kinetic 
energy T att under the left action of Q that it can be expressed as a (positive, symmetric) 
bilinear form in p and the left-invariant angular velocities Ui and Q a (i,a = 1,2,3). The 
bilinear form is further restricted by the identifications ([4.11|) which imply that we should 
identify the the left-invariant angular velocities uj and ft calculated from D(C) and D(G) 
with those calculated from the image of D(C) and D(G) under the right action of On, 0\ 2 
and Oq3- Explicitly, the angular velocities transform as follows: 



0.3 



O 
On 

Ol2 



(fii, Q 2 , ^3) i-> (^1, ^2, ^3) and (uji, u 2 , uj 3 ) i-> (-a>i, -u 2 , uj 3 ) 

(fii, tt 2 , fi 3 ) 1 ► (fii, -0 2 , -^3) and (wj, cj 2 , cj 3 ) i-> (u u -u 2 , -u 3 ) 
(fii, f2 2 , fi 3 ) 1 (fii, -0 2 , -fi 3 ) and (wj, cj 2 , cj 3 ) i-> (-wj, cj 2 , -o; 3 ). 



The kinetic energy must respect these identifications, which implies that it is of the following 
form in the centre of mass frame (S = 0): 

T M = ^ [fp 2 + a 2 u)\ + b 2 u 2 2 + d 2 u 2 3 + A 2 tt\ + B 2 Q 2 2 + C 2 Q 2 3 + 2eu 3 Q 3 ) , (5.3) 
where a, b, d, e, /, A, B and C are functions of p only. Explicitly we find 

f 2 = J d 3 x {-tv(R p R p ) - ^tr([i^][£ p ,4])} • (5.4) 

For the spatial moments of inertia we have 

a 2 = J d 3 x |-tr(aj x R)\ - ^tr([(aj x 4] [(x x A)^])} (5.5) 

and similar formulae for 6 2 and d 2 with the free index 1 replaced by 2 and 3 respectively. 
The moments of inertia in iso-space are given by 

A 2 = J rf 3 x{^tr([r 1 ,f>]f>t[r 1 ,f>]m) + ^tr([[r 1 ,f>]m,i? l ][[r 1 ,^m,4])} (5.6) 
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and similar formulae for B 2 and C 2 with t\ replaced by r 2 and t 3 respectively. Finally the 
coefficient e of the cross term is 

e = ijd 3 x {-±tr([ T3 ,ft]fr(x x R) 3 ) - itr([[r 3 , U]U\ Ri) [{x x R) 3 ,Ri))} . (5.7) 

In all the above integrals the integrands are invariant under the maps Oq 3 ,Oh,Oi2- As a 
result they can be computed by integrating only over the octant X{ > (i — 1, 2, 3) and 
multiplying the result by 8. Details of how this is done in practice may be found in the 
appendix. 

When S / the only additional term in the kinetic energy that is allowed by the 
requirement of invariance under O 03 , On and Oi 2 is 

DWikDWflSiSjMu (5.8) 
where the tensor M^i depends on p and can be calculated from 

M M = J d 3 x |-itr(itR fc ) " ~tr([A z A][4A])} • (5-9) 

The reflection symmetries of Ri imply that Mik is diagonal, so that one needs to calculate 
three more functions of p in order to study the dynamics in an arbitrary inertial frame. 
In the following we will restrict attention to the centre of mass frame. Thus we will only 
compute the Lagrangian on the manifold M7 of centred attractive channel fields. 

It is useful to rewrite the expression ( p.3| ) slightly by completing the square: 

T att = ^ (f 2 p 2 + a 2 u\ + bW 2 + c 2 uj\ + A 2 n\ + B 2 Vl\ + C 2 (Q 3 + wu 3 ) 2 ) , (5.10) 

with c 2 = d 2 — e 2 /C 2 and w = e/C 2 . The kinetic energy can be expressed in terms of the 
metric 

g = fdp 2 + a 2 o\ + b 2 o\ + c 2 o\ + A 2 Y,\ + B 2 Y? 2 + C 2 (S 3 + wa 3 ) 2 (5.11) 

where <7j and S a (i, a = 1,2, 3) are the left-invariant one-forms which, when evaluated on a 
tangent vector to a trajectory in M 7 , give the left-invariant angular velocities: 

a i(ll) = Ui and = ^o- ( 5 - 12 ) 

dt dt 

Explicitly, one finds the following formulae in terms of the Euler angles defined earlier 

o\ = sin ipdO — cos ij) sin Odcf) 
o"2 = cos ipdd + sin ip sin Odcf) 

a 3 = dip + cos 9d<f), (5.13) 
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and the corresponding expressions for S a with (<f>,9,ip) replaced by (a, (3, 7). One checks 
that doi = —\tijk i A (Tfc and dT, a = — |e a fe c S fe A S c . For our discussion we also require the 
left- invariant vector fields £j on 50(3) J and on 50(3) 7 which are dual to the forms a,i 
and S , i.e. = Sy and £ a (G>) — <W Explicitly, in terms of Euler angles, 

8 8 cos ip 8 

4i = cot 6^ cos 777 + sin ■0- 



dip 89 sin# d(f) 

8 .8 sin ip 8 

62 = — cotOsmip— + cosip— + 



dip 86 sin 9 8<j) 

& - I- < 514 > 

The corresponding formulae for the vector fields (b are obtained by replacing (0, 6, ip) with 

(a, (3, 7). 

For general p, the determination of the metric g requires the computation of eight func- 
tions of p, but when p = p the form of the metric is constrained by the invariance of the 
field U(p ,x) under the 5*0(2) right action ( |4.8|) . It follows that we should identify tangent 



vectors at po which are related by the 5*0(2) right action ( |4.12| ). That action can be used 
to rotate £1 into £2 and £i into £2- Since the kinetic energy must respect this identification 
we conclude that, for p = p , #(61,61) = #(6,62) and #((1,(2) = #((2,(2), which implies 

a 2 (p ) = b 2 (p ) and A 2 (p ) = B 2 (p ). (5.15) 

Furthermore, the vector field 63 + 2C3 generating the 50(2) action is not defined when p = p , 
and for the metric to be regular there we require that 

g(& + 2(3,6 + 2Cs) = c 2 (p ) + 2 ( Po )(2 + w(p )) 2 = 0. (5.16) 

Thus, assuming O(p ) 7^ (equality would imply that ( 3 also has length zero), it follows 
that 

c(p ) = and w{p ) = -2. (5.17) 

For large p, U(p, x) is approximately of the product form, and the potential and kinetic 
energy can be expressed in terms of the mass and the moment of inertia of the hedgehog 
field Uh (|3.5|) and the parameter A characterising its profile function. Such a calculation 
was carried out in for the product ansatz constructed from Lorentz-boosted hedgehog 
solutions. We have repeated that calculation using the instanton-generated hedgehog field for 
the individual Skyrmions and omitting relativistic corrections. Both of these modifications 
actually make the calculations harder: in [|15| one could exploit the fact that the hedgehog 
fields in the product ansatz individually satisfy the static equations of motion, and it was 
found that the relativistic corrections make the final answer simpler than it would otherwise 
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be. Here we take the centred product ansatz for the attractive channel ( 3.10|) and act with 
spatial and iso-spatial rotations to obtain 

CU H {D(G)- 1 x + iRe^UniDiG)-^ - iRe 1 )r 3 C jf . (5.18) 

Then we insert this field into the Skyrme Lagrangian, allowing the SU(2) matrices C and G 
and the separation parameter R to vary with time. 

The results can be expressed in terms of a single Skyrmion's mass, M = 147.2, its moment 
of inertia, A = 140.1, and its dipole strength, p = 3.31. For the potential energy we find the 
asymptotic formula 

V~2M-1M — (5.19) 
R 6 

and for the moments of inertia 

a 2 ~ 2A 

b 2 ~ lMi? 2 + 2A-1.2-^ 

d 2 ~ lMi? 2 + 2A-0.16-^ 
1 R 

A 2 ~ 2A 

47Tp 2 



B 2 ~ 2A + 2.0- 



C 2 ~ 2A-2.0- 



i2 

47Tp 2 

I 



47T» 2 

-2A + 0.64-— (5.20) 



with corrections of order 1/R 2 . It follows that 

z it 
w 1_ .68-^, (5.21) 

also with corrections of order 1/R 2 . Apart from the 1/R corrections the formulae for the 
moments of inertia can be understood quite easily in terms of the moments of inertia of 
the individual Skyrmions and Steiner's theorem. Finally the asymptotic form of the radial 
metric coefficient is 

/ 2 (p)~iM + 2.72-^. (5.22) 

As mentioned in section 3, the parameter R can be identified with p for large R. Thus the 
above formulae can be compared with the asymptotic form of the numerically calculated mo- 
ments of inertia. The agreement is very good for p > 10 . In figures 2 - 6 we plot cubic splines 
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constructed from our numerical values for the metric coefficients a 2 , b 2 , c 2 , A 2 , B 2 , C 2 , w and 
f 2 ; the values for p > 15 were calculated using the asymptotic expressions ( |5.20p - ( [5.22j ). 



The numerical values at p = p are a 2 = 348.5, b 2 = 348.4, c 2 = 0.0, A 2 = B 2 = 229.8, 
C 2 = 130.3, w = -2.0 and f 2 = 199.9. 

6 Quantisation 

Our next goal is to write down the quantum Hamiltonian in the attractive channel and in 
the centre of mass frame. As a first step we calculate the classical Hamiltonian from the 
Lagrangian L att = T att — V in the usual way. Thus we define the momenta conjugate to the 
angular velocities u> and O, 

which are the body-fixed angular momenta in space and iso-space. In components, 

L 1 = a 2 uj 1 , L 2 = b 2 uo 2 , L 3 = (c 2 + C 2 w 2 )u 3 + wC 2 £l 3 , , 

K 1 =A 2 Q 1 , K 2 = B 2 Vl 2 , K 3 = C 2 (n 3 + wuo 3 ). lb ] 

The corresponding space-fixed spatial angular momentum is J = D(G)L and the space-fixed 
iso-spatial angular momentum is / = D(C)K. Finally defining P p = dL M /dp we arrive at 
the classical Hamiltonian in the centre-of-mass frame 

The conserved quantities are the Hamiltonian itself, J and /. It follows that L 2 = J 2 and 
K 2 = I 2 are also conserved. 

In our quantisation scheme the quantum Hamiltonian is 

h 2 

H = — —A + V. (6.4) 



Here A is the covariant Laplace operator associated with the metric g ( |5.11| ). Explicitly, in 
terms of the vector fields and Q defined in the previous section, 

A = 1 9 ( abcABC d \ g g - ^C 3 ) 2 Cf , Cf , Q ( , ~ 
abcABCf dp\ f dp) a 2 P c 2 A 2 B 2 C 2 ' 1 J 

Physically, —ih^i and —ihQ are the components of the operators for the angular momenta 
L and K respectively. 

For the rest of this paper we will be concerned with the stationary Schrodinger equation 

= EV. (6.6) 
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The wavefunction \I/ is a section of a (possibly trivial) - complex line bundle over M7. The 
structure of the bundle is fixed by Finkelstein-Rubinstein constraints, which we summarise 
as follows (see also the discussion in Q). If a single Skyrmion is quantised as a spin | particle 
then the wavefunction for several Skyrmions should pick up a minus sign when one of the 
Skyrmions is rotated by 2ir. Since a rotation of a single Skyrmion by 2it is a non-contractible 
loop in Qb for any B G Z, and since moreover all non-contractible loops are homotopic to 
this one, wavefunctions are required to be sections of a line bundle over Qb whose holonomy 
around any non-contractible loop in Qb is — 1. It follows from the fact that the exchange of 
two Skyrmions is a non-contractible loop that Skyrmions quantised as half-odd-integer spin 
particles are fermions and Skyrmions quantised as integer spin particles are bosons. This 
result, due to Finkelstein and Rubinstein ||, is an example of a topological spin-statistics 
theorem. 

In the section 4 we showed that M7 has the same zeroth and first homotopy groups as 
Q 2 - Thus, since we quantise a single Skyrmion as a spin ~ particle, ^ should be a section 
of a non-trivial bundle over M 7 such that its holonomy around a non-contractible loop in 
M7 is —1. It is sufficient to impose this for one such loop, and we will use the loop r n , 
which, for well-separated Skyrmions, is the rotation of one Skyrmion by 2n. The simplest 
way to implement the Finkelstein-Rubinstein constraints is to think of $ as a function on 
the double cover M7 of M7 and to impose the equivariance condition 

*oOu = -*, (6.7) 
Points related by O03 are still identified in M 7 , so \I/ must obey 

$oO 3 = *. (6.8) 
Then, since O i2 = O11O03 it follows that 

$oOi 2 = -^. (6.9) 

In accordance with our interpretation of 0\i this last condition shows that the two Skyrmions 
are fermions, i.e. that the wavefunction is odd under their exchange. The important point 
here is that M 7 captures enough of the topology of Q2 for the topological spin-statistics 
theorem to hold. 

Exploiting the symmetry of the Hamiltonian we separate variables by introducing Wigner 
functions D J sm ((j) , 9 , ip) and -D£ n (ai, /3, 7) on SO(3) J and SOft) 1 respectively, following the 
conventions of [16]. The former are simultaneous eigenfunctions of the square of the total spin 
operator L 2 = J 2 = —h 2 (^f + £f + £f ), the third component of the space-fixed spin operator 
J 3 = —ihd/dcj) and the third component of the body-fixed spin operator L 3 = —ifid / difj: 

L 2 Dl m = j(j + l)h 2 Di m , J 3 Dl m = mhDi m , UD{ m = shDi m . (6.10) 
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Similarly, D\ n (a, /3,j) is an eigenfunction of the square of the total iso-spin operator K 2 = 
I 2 = — ^ 2 (Ci+C2+Cf)> the third component of the space-fixed iso-spin operator ^3 = —ihd/da 
and the third component of the body-fixed iso-spin operator K% = —ihd/d^f. 

K 2 D\ n = i(i + l)hD\ n , hD\ n = nhDl, K,D\ n = thD\ n . (6.11) 

The Wigner functions are normalised so that 

/ \D{ m \ 2 sin 9 d(f> dtp d9 = and f \D l tn \ 2 sin (3 dad-? dp = ^—r- (6.12) 

J ^7 I 1 J I 1 

Since both spatial and iso-spatial rotations by 2tt are contractible in M 7 we must use 
Wigner functions which are even under such rotations and hence i and j are positive integers 
and n, t, m, s are integers in the ranges —i<n,t<i and — j < m, s < j . Wigner functions 
of integer spin transform under the discrete maps ( |6.8| ) and ( |6.7| ) as follows 



DL°O 03 = (-l) s Dl m 
D\ n Di m oO n = { 1 )• •'/;' ,,./>',,,„. (6.13) 

When looking for bound states of the Schrodinger equation ( |6.6|) we can fix the quantum 
numbers j (total spin), i (total iso-spin), m (the third component of the space- fixed spin) 
and n (the third component of the space-fixed iso-spin). We are interested in two nucleon 
bound states, so i = or i = 1. The Wigner function for i = j = (which is constant) does 
not satisfy ( |6.13|) . For the quantum numbers of the deuteron, = (0, 1), there is exactly 



one state (we do not count Wigner functions differing in m or n as different states), which 
we write as 



* = y^L(0,0,vo«(p)- (6-i4) 

Inserting \l/ into ( |6.6|) we obtain the following differential equation for u: 

1 d ( abcABC du\ (1, 1 1\ . 

X + l^(E-V)-- 2 -- )« = 0. (6.15) 



abcABCfdp\ f dp) Vr v ' a 2 b 2 . 

For = (1, 0), the iso-vector states, there are two allowed angular states. The ansatz 



* = J-^Dl n (a,l3,i) U (j>) (6.16) 



57T 2 " 



leads to the differential equation 



1 d ( abcABC du\ 2 , > 1 1 , 

+ (^(E-V)---- « = 0. (6.17) 



abcABC f dp \ f dp \h 2K A 2 B 2 
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The second possibility is 



and leads to 



* = V + D \n)(a, P, l)u{p) (6.18) 



d (abcABCdu\ ( 2 . „ TrX 1 w 2 1 \ 



abcABCfdp\ f dp J \h 2 " 'A 2 c 2 C 2 



To study more general bound states and scattering solutions of ( |6.6|) one can use tech- 
niques similar to the ones used in |1| to study the quantum dynamics of two magnetic 



monopoles in the moduli space approximation. One should make the ansatz 

tt(0,M,a,/8,7,/>)= E £>L(a, P, 7 )£>L(0, 0, V0<Up), (6-20) 

where the sum runs over the indices of the Wigner functions, suitably restricted by conser- 
vation laws and the constraints (|6.7| ) - ( |6.9| ). For bound states, one may fix i, j, n and m, 
and one will typically need to solve a set of coupled ordinary differential equations for the 
functions Ut nsm (p), —i<t<i,—j<s<j. To obtain scattering solutions one necessarily 
needs to consider an infinite sum over j, and in principle one needs to solve infinitely many 
systems of coupled ordinary differential equations. However, to compute scattering cross 
sections at low energy it is sufficient to study the systems of equations that arise for small j. 

In this paper we will only study bound state problems and restrict attention to the 
uncoupled ordinary differential equations ( |6.15|) , ( |6.17|) and (|6.19| ). 



7 Quantum Bound States 



To find the lowest eigenvalue of the deuteron equation ( |6.15| ) we need approximate analytic 



solutions near p = po and for large p. To find these it is useful to define the effective potential 

M _ 2M + f (1 + I-^L), (7.1) 

which tends to zero as p — ► oo. In figure 7 we plot both the classical potential V — 2M and 
the effective potential V eB . Also defining the shifted energy 

e = E-2M-^, (7.2) 



equation ( |6.15| ) becomes 



1 d ( abcABC du\ 2 , Tr , . . 

+ -^{e-V ce )u = Q. (7.3) 



abcABCf dp \ f dp J h 
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For p near p , we define h = p — po , and the equation becomes approximately the Bessel 
equation 



hdh\dh) h 



h ~ ) + -^tF^ (e - K a (po)) « = 0, (7.4) 



which has only one solution that is regular at the origin, namely the Bessel function Jo of 
zeroth order. Thus, for small h 



u( Po +h)ocJ [ tp^^ (e - K ff (po))/i ] • (7.5) 



It follows in particular that u is non-zero at po and that it has vanishing derivative there. 
For large p the equation ( |7.3|) becomes 

ld 2 (pu) Me n 

which has two solutions, one exponentially growing and one exponentially decaying. A 
normalisable wavefunction has to be asymptotically proportional to the latter: 

e — p 

u oc . (7.7) 

P 

Imposing the behaviour ( [7.5|) and ( [7.7|) near p and for large p we have numerically searched 
for normalisable eigenfunctions of ( [7.3|) , using a shooting method; we find that there is 
precisely one, which we denote by Ud and which we normalise so that 

dpu 2 d abcABCf = 1. (7.8) 

PO 

The corresponding eigenvalue is Ed = —1.107 in geometrical units, or Ed = —6.18 MeV in 
physical units. We also note that the corresponding eigenvalue of the original Schrodinger 
equation (gjp is, according to ( ffT2D , E d = 2M + h 2 /AA + E d = 1658.4 MeV in physical units. 



In the following, the wavefunction ( 5.14 ) with Ud for u and m—1 ('spin up'), namely 



^ = y^oi(<M^K(p) = -^^^Be**u A {p) (7.9) 

will be referred to as the deuteron wavefunction. In figures 8. a) and 8.b) we show plots of 
Ud and of the normalised probability density u 2 d (p)abcABC f . 

The translation of the eigenvalue Ed of ( |7.3|) into a theoretical prediction for the deuteron 
binding energy requires some thought. Physically, the deuteron's binding energy is defined 
as the energy required to break a deuteron up into two infinitely separated nucleons, one 
proton and one neutron. However, in this paper we are only considering two-Skyrmions in the 
attractive channel, and therefore wavefunctions on M-?, the space of attractive channel fields, 
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can only describe two nucleons with correlated spin and iso-spin. Thus, while the energy of 
two free nucleons is 2M + 3fr 2 /4A = 1708.2 MeV, the eigenvalue e d measures the binding 
relative to the energy of infinitely separated two nucleons in the attractive channel and with 
total iso-spin and spin quantum numbers = (0, 1), which is 2M + fr 2 /4A or 1664.6 MeV. 
However, it would be naive to simply add the difference of —h 2 /2A = — 43.6MeV to e d and 
interpret the result as the true binding energy of the deuteron. Rather, a precise calculation 
of the deuteron's binding energy would require solving the Schrodinger equation on the full 
twelve-dimensional moduli space Myi discussed in the introduction. That is a difficult task, 
but it may be possible to take into account the extra degrees of freedom approximately by 
treating them as small oscillations around the attractive channel fields. Then one might be 
able to estimate the effect of the extra degrees of freedom on the Schrodinger equation for 
the deuteron by including their zero point energies in the effective potential V cB . These zero 
point energies depend on p, but in a first, crude approximation one may describe the net 
effect by adding a constant to V cS . From the above discussion of well-separated nucleons in 
our model we know that the constant should be approximately 43.6 MeV. It then follows that 
the binding energy of —6.18 MeV, which we calculated by restricting attention to attractive 
channel fields, is also an estimate of the binding energy calculated in a more careful treatment 
involving the larger manifold of collective coordinates Myi- 

It is interesting to convert the quantum mechanical probability distribution |\I/| 2 (p, 0, 9, i/j) 
of a quantum state \l/ on the moduli space into a probability distribution pq,(x) on physical 
space in such a way that pq,(x)d 3 x can be interpreted as the probability of finding one of the 
nucleons in the region [xi,xi + dxi] x [x 2 , x 2 + dx 2 ] x [x 3 , x 3 + dx 3 ]. Strictly speaking such a 
translation only makes sense for Skyrme fields consisting of two well-separated Skyrmions. 
For such fields the baryon number density B° is peaked near the individual Skyrmions' 
positions and may be interpreted as the 'particle density' of the Skyrme field. When the 
Skyrmions coalesce it is no longer meaningful to talk about the individual Skyrmions, but 
one can still interpret the baryon number density as the 'density of Skyrme matter'. To 
find the probability distribution p^(x) of Skyrme matter for a given quantum state \I> on 
the moduli space, one should average the classical baryon number density over the moduli 
space, weighted by \^\ 2 - Thus, writing B°(p,x) for the baryon number density ( |2.12j ) of 
an attractive channel field in the standard orientation and with separation parameter p, we 
define the spatial probability distribution Pd{x) for the deuteron state as 

p d (x) = -j B°{ P ,D(G)- 1 x)\^ d \ 2 sin 9 abcABC f dp d0d<f> #. (7.10) 

Here the SO (3) matrix D(G) is parametrised in terms of the Euler angles (0, 9, ip) according 
to (|4.13|) . To evaluate this expression we introduce spherical coordinates (r, 6, $) for x, so 
that x = r(sin B cos $, sin B sin $, cos 6) and a further set of spherical coordinates (r,9,<j>) 
for x = D(G)~ 1 x, i.e. x = r(sin6 l cos0, sin 9 sin 0, cos#). Then, by expanding B°(p,x) for 
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fixed p and r in spherical harmonics Yi m (9, 0), and by using the transformation properties of 
spherical harmonics under rotations, one shows that p d depends only on r and 0, and that 
it is given by 

p d (r,e) = —f B°(p,x)u d (p)abcABCf sin 6 dp Sdi 

87T J 

— (1-3 cos 2 6) / B°(p,x)u 2 d (p)(l -3 cos 2 6) abcABCf sin 6 dp d9d4>. 
32% J 

(7.11) 



In the conventional description of the deuteron as a bound state of two point-like nucleons, 
the square of the modulus of the deuteron wavefunction also gives the probability of finding 
one of the nucleons in a given region of space. It is therefore meaningful to compare pd{r, 0) 
with such a probability distribution calculated in a conventional model. That comparison 
is made in figure 9. The distributions are clearly very similar. Thus, although the classical 
toroidal bound state of two Skyrmions looks radically different from the conventional picture 
of the deuteron (a point often raised as a criticism of the Skyrme model) the two approaches 
lead to remarkably similar spatial probability distributions of nuclear matter at the quantum 
level. 

To study the quantum states with = (1,0) we need to consider the two radial 

equations (|6.17| ) and ( p,19| ). There is no bound state solution of ( |6.19|) , essentially because 



c 2 vanishes at po, giving rise to a strongly repulsive centrifugal potential w 2 /c 2 . The equation 
(|6.17 ), however, does have a bound state solution. It can be brought into the same form as 



the deuteron equation (|7.3j), but now the shifted energy e is 



2 



E-2M- — (7.12) 
2A v ' 



h 2 ( l l i 



and the effective potential is 

v -= v - 2M+ 'i{^ + w-i)^ <" 3 > 

which tends to zero as p — > oo. There is a unique bound state solution which, like the 
deuteron wavefunction, is non-vanishing at po and has zero derivative there. We will not go 
into the details here, but simply note the numerical value of the bound state energy, which 
is Enn = —1.74 in geometrical units, or enn = —9.74 MeV in physical units. It follows that 
E NN = 1676.7 MeV. 

In nature, the deuteron is the only bound state of two nucleons. However, there is a 
marginally unbound state with quantum numbers = (1,0), the iso-vector 1 S'o state. In 
our model this state is the one with energy Enn- Our model gives the right ordering of 
energy levels E d < E NN , and an energy difference E NN — E d of the right order of magnitude. 
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Moreover, Enn + 43.6 MeV > 1708.2 MeV, the energy of two free nucleons. Thus if, as 
argued above, the effect of including the two extra degrees of freedom is to push the energy 
levels up by 43.6 MeV, then the = (1,0) bound state we found here will not persist 
when all 12 degrees of freedom of M 12 are taken into account. 



8 Electrostatic Properties of the Deuteron 



The first attempt to calculate the Skyrme model's predictions for the deuteron's electrostatic 
properties was made by Braaten and Carson in under the assumption that the deuteron 
can be described as a quantum state of a toroidal Skyrmion. We will briefly review and then 
use their basic formulae for the classical electrostatic properties of a Skyrme field. Since we 
are using a larger number of collective coordinates than in M, our calculation of quantum 
mechanical expectation values will, however, be different. 



The starting point of the calculations is the basic relationship, first derived in [£L7|], be- 



tween the electromagnetic current j M , the third component of the (classical) iso-spin current 



and the baryon number current introduced in ( 2.11|) : 



3 ^ = \B, + Il. (8.1) 

We are only interested in the expectation values of various classical multipole tensors for 
quantum states with zero iso-spin, so we set to zero in the following and replace j M by 
\B lL . It is sufficient to perform the computations of the various tensors with the Skyrme 
fields in their standard orientation. Specifically, we will compute the theoretical predictions 
for the deuteron's electric charge radius, its electric quadrupole moment and its magnetic 
dipole moment. We will first discuss the general formulae for these quantities and then 
compile the numerical results in a table. 

First consider the root mean square (rms) electric charge radius of an attractive channel 
field. It is a function of p only and defined as the square root of 



r 2 

rms 



(p) = lJd 3 x\x\ 2 B°(p,x). (8.2) 

The theoretical prediction for the deuteron's rms electric charge radius r c is the square root 
of 

POO 

rl = <* d |rLl* d > = / dprlu 2 d abcABCf. (8.3) 

Next, the classical electric quadrupole tensor for an attractive channel field in standard 
orientation also depends only on p and is given by 

Qij(p) = \ J (Pxi&CiXj - \x\%)B°{p, x). (8.4) 
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It follows from the discrete symmetries of the attractive channel fields that Qij is diagonal. It 
is also traceless, so it contains only 2 independent components (Qu and Q 22 , say). Moreover, 
at po, Qu = Q 22 - The quadrupole tensor for a field in general orientation G is given by 

Q lm = D(G)uD(G) mj Qij (8.5) 

and the deuteron's quadrupole moment is defined as 

Q=(*d|Q33|*d)- (8.6) 
Performing the angular integration one finds 

Q = \ dp (Q n + Q 22 )u 2 d abcABCf. (8.7) 

Finally, the iso-scalar (I 3 = 0) part of the magnetic dipole moment is defined as 

fJ'i(p) = \ J e ijk x j B k (p,x)d 3 x. (8.8) 
For fields in the standard orientation it can be written 

Pi = M ia VL a + rh ik uj k , (8.9) 

where 



m ik = --^ / d 3 xe kmn x j x m tT([R i , RjjRn). (8.10) 

It follows from the discrete symmetries of the attractive channel fields that both Mi a and 
rhu are symmetric, and that the only non-zero component of M{ a is M33. In fact, we shall 
see shortly that we only need rhu and rh 2 2 when computing the deuteron's magnetic dipole 
moment. Explicitly one finds 

rhn(p) = \ J d 3 %(4 + 4)B°(p, x) 

m 22 (p) = \ J d 3 x(xl + x 2 2 )B°(p, x); (8.11) 

at po these two quantities are equal. To compute the deuteron's magnetic moment, defined 
via 

p= (* d \ti 3 \* d ), (8.12) 

we require the magnetic moments of the attractive channel fields in arbitrary orientation G, 
which are given by 

p^p) = D(G) H p t (p). (8.13) 
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Then, in order to calculate the expectation value in the quantum state of the deuteron, we 
should replace ct>i,a>2 and 0J3 by the operators Li/a 2 , L 2 /b 2 and L 3 /c 2 , and similarly fii, Q 2 
and ^3 by the operators Ki/A 2 , K 2 /B 2 and K 3 /C 2 . Since the components of the body-fixed 
angular momentum operator L do not commute with the entries of the matrix D(G) there 
is potentially an operator ordering ambiguity in evaluating matrix elements of the operator 
fjLi, but in the calculation of the expectation value in the deuteron state this ambiguity does 
not arise. Here we find 

^=U C °dp(^f + ^)u 2 d abcABCf. (8.14) 

A J pq (X 

In nature, the magnetic dipole moment of the deuteron is almost exactly equal to the sum 
of the magnetic dipole moments \i p and \x n of the proton and the neutron. In conventional 
models of the deuteron as a bound state of a proton and neutron the discrepancy can 
be related to the d-wave contribution in the deuteron wavefunction It is instructive 



for us to carry out a similar comparison, but here we should compare our result for the 
deuteron's magnetic dipole moment with the sum of the proton's and neutron's magnetic 
dipole moments as calculated in the Skyrme model with the appropriate parameters. 

We list our numerical results in the table below. We have followed the practice in nu- 
clear physics of measuring magnetic dipole moments in units of the nuclear magneton (nm) 
h/2M^, where = 939 MeV is the physical nucleon mass. For comparison we also list 
the results obtained by Braaten and Carson in their treatment of the deuteron as a toroidal 
Skyrmion. The experimental values are taken from [ZHfl. In the first row we compare val- 



ues for the deuteron binding energy. Under Braaten & Carson we quote the value which 
Braaten and Carson calculated in Q by computing the difference between the energy of 
their deuteron state and the energy of two free nucleons. Braaten and Carson did not in fact 
attach much meaning to this particular result of their calculations and pointed out that one 
should take into account vibrational modes of the toroidal configurations for a meaningful 
calculation of the deuteron's binding energy. The comparison in our table shows that the 
inclusion of the softest vibrational mode alone - the deformation of a toroidal configuration 
within the attractive channel manifold M 10 - reduces the prediction for the binding energy 
by an order of magnitude, and brings it close to the experimental value. In the last row, the 
theoretical predictions for \i v and \i n are calculated in the Skyrme model using the formulae 
given by Adkins et al. in ||. These formulae depend on the profile function for the hedgehog 
field with unit baryon number and on the choice of the constants F n and e. For the latter 
we use the values assumed throughout this paper and given in (|2~5| ). Then, to calculate the 
value of Hp + /i n in our model we use the instanton generated profile function ( |3.6|) ; the value 
listed under Braaten & Carson is taken from the analysis of Adkins and Nappi in J/J which 
is based on the same values for F n , e and the pion mass as used by Braaten and Carson. 
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Theory 


Experiment 




this model 


Braaten & Carson 




e d [MeV] 


-6.18 


-158 


-2.225 


r c [fermi] 


2.18 


0.92 


2.095 


Q [fermi 2 ] 


0.83 


0.082 


0.2859 


jj, [nm] 


0.55 


0.74 


0.8574 


/i p + fi n [nm] 


0.41 


0.55 (from §) 


0.8797 



Table 1 

Deuteron Properties 



9 Discussion 

In this paper we have described the deuteron as a quantum state of instanton- generated 
two-Skyrmions in the attractive channel. This picture seems qualitatively correct and leads 
to predictions for certain deuteron observables which are also in reasonable quantitative 
agreement with experiment. However, our calculations clearly involve a number of approx- 
imations, and from the point of view of the general approach to two-Skyrmion dynamics 
outlined in the introduction, two of these deserve a more careful discussion. Firstly, we 
obtained our Skyrme fields by computing instanton holonomies and not, as suggested in the 
introduction, by calculating paths of steepest descent in the Skyrme model. This means in 
particular that we cannot take into account the physical pion mass. The experience with 
the static properties of a single nucleon in the Skyrme model shows that the inclusion of the 
pion mass does not affect the results very much. We expect similarly that the plots 2-6 and 
7 of the metric coefficients and the potential would not be qualitatively different for fields 
calculated via paths of steepest descent in the Skyrme model with a pion mass term. There 
is one obvious difference which would occur in the asymptotic form of the plots: while the 
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potential and metric coefficients calculated above approach their asymptotic values accord- 
ing to a power law, the approach would be exponentially fast in a model with non-zero pion 
mass. This difference would be important if one were to discuss nucleon-nucleon scattering. 
However, it does not matter much for the calculations presented above because the deuteron 
wavefunction is not very sensitive to the asymptotic form of the Hamiltonian. 

The other, more serious simplification we made was to consider the ten-dimensional 
moduli space Mio instead of the twelve-dimensional space M\2 described in the introduction. 
The validity of this approximation is much harder to assess, but the following comments may 
elucidate its physical meaning. It has been known for a long time (see e.g. that, for 

zero pion mass, the potential energy for the interaction of two well-separated Skyrmions with 
arbitrary individual positions and orientations is, up to a positive constant of proportionality, 

A-3(n-R) 2 . . 

-(i-cosx) ! ^- J -, (9.1) 

where R is the relative position vector, R = R/R and (n,x) are the axis and rotation 
angle for the SO (3) matrix describing the relative orientation. In the attractive channel, the 
relative orientation is chosen so as to minimise the potential energy at fixed R. Thus, \ is 
set to 7r and n is required to be orthogonal to R. This eliminates two degrees of freedom, 
but it still allows for rotations of n in the plane orthogonal to R. 

The potential ( |9.1| ) is the classical analogue, in the Skyrme model, of the one pion ex- 
change tensor potential in conventional nuclear physics ||22|| . To calculate nuclear forces from 
it one should compute its expectation value in states which are tensor products of the free 
nucleon wave functions given in 0. In particular, the expectation value of ( |9.1| ) in the triplet 
state with total spin 1 and total isospin (the quantum numbers of the deuteron) is a 3 x 3 
matrix which can be expressed in terms of the total spin operator s. It is, up to an overall 
positive constant of proportionality 

V^ 2 -^- (9.2) 

At fixed R, this expression is minimal when the total spin is parallel to the separation vector 
R. The assumption of the attractive channel approximation is that the relative orientation is 
always such that the potential is minimal at given R. Thus, working with attractive channel 
fields amounts, in the language of nuclear physics, to assuming that the total nucleon spin is 
always aligned with the relative separation vector or, equivalently, that the torque resulting 
from the tensor force is infinitely strong. 

In conventional discussions of the deuteron and its properties the tensor force is respon- 
sible for the existence of a d-wave contribution to the deuteron wavefunction The d-wave 
probability in turn is linked to physical observables. In the absence of a d-wave the deuteron's 
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electric quadrupole moment would be zero and the deuteron's magnetic dipole moment would 
equal the sum of the proton's and the neutron's magnetic dipole moments. Both the size of 
the quadrupole moment Q and the difference /i — + /i n ) are therefore direct measures of 
the d-wave probability, which in turn indicates the strength of the tensor potential. Thus, 
since the truncation of M 12 to the space of attractive channel fields systematically overesti- 
mates the strength of the tensor force, it is not surprising that our theoretical predictions for 
Q and \i — (p p + p n ) are rather large. Quantising the extra two degrees of freedom included in 
M12, at least approximately, may well bring the predictions much closer to the experimental 
values. 
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A Numerical Methods 

In this appendix we give a detailed and in parts technical account of our method for comput- 
ing the potential V in the attractive channel, the metric g derived from the kinetic energy 
T att , and the quantities required in the calculation of static deuteron properties. All these de- 
pend only on the Skyrmion separation p, and in equations (|2.10|) , ( |5.4| )- ([5T7| ) and (|3.2|) - (|8.11 



we have shown how to write them as integrals over physical space of various combinations 
of the Skyrme fields U (p, x) and their currents. The metric g has eight non- vanishing com- 
ponents, and it turns out that the deuteron's static properties can be expressed in terms of 
three independent moments of the baryon number density B°. Including the potential, there 
are thus twelve functions of p to be computed. 

The computational method is similar to that used in [E3| to find classical bound states of 



three and four Skyrmions, although several additions and refinements are necessary to carry 
the programme through successfully. The following paragraphs concentrate on these new 
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features; [^3] should be consulted for the remaining details. There are two tasks: firstly the 
solution of the holonomy equation ( |3.2| ) to construct the attractive channel Skyrme fields in 
standard orientation, and then integrations over all space to obtain the required functions 
of p. 

Recall that the attractive channel Skyrme fields are obtained from instantons which 
depend on two essential parameters. In the standard orientation, for which the relevant 
Hartshorne data is depicted in figure 1, these may be taken to be the scale L and the angle 
d. Thus, the solutions of the holonomy equation (|3.2| ) for these instantons are Skyrme fields 
on R 3 which also depend on L and we denote them by U(L, i?, x). It follows from the 
definition of U (L, i?, x) and from the scaling behaviour of JNR instantons that 

U(L,#,x) = U{l,#,j). (A.l) 

Hence the current Ri(L, x) = (diU(L,'d,x))U'(L,'d,x) satisfies 

R i (L,#,x) = jR-(l,#,j), (A.2) 

where on the right-hand side the superscript L indicates that the differentiation should be 
carried out with respect to x/L. To obtain the Skyrme fields U(p, x) used in the main body 
of the paper, the scale L is fixed for each value of $ at the value L(d) which minimises 
the potential energy ( [2.10D . Explicitly, writing E 2 and E± for the quadratic and the quartic 
terms in the potential energy ( |2.10|) evaluated on the field £/(l,$, x), this is 

m = M. (a.3) 

V -^2 

In practice, it is most convenient to compute all fields and currents at the scale L = 1 and 
to find the fields and currents at the relevant scale using the formulae ( A.lj) and ( A.2 ). 

It remains to explain the computation of the current R p (|5.1|) , which requires some care. 
After fixing the scale at Li^d) the formula p = 2L(i9)(l — sin 1 *?) establishes a one-to-one 
relation between p and Thus, both •& and L may be thought of as functions of p, and we 
can write 

d p U(L, 0, x) = (^)d L U(L, 0, x) + (^)a^(L, x). (A.4) 

Using the relation ( |A.l ), the derivative of U with respect to L may be expressed in terms of 
the currents Ri\ 

d L U(L, x)U\L, &, x) = - E(t) x )- ( A - 5 ) 

i L 

Thus it is the derivative of U with respect to or equivalently the current R$ = (d#U)W, 
that has to be calculated directly from the holonomy equations. Finally, the linear combi- 
nation on the right-hand side of ( |A.4j ) is computed after the scaling factors £(t?) have been 
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calculated; dL/dp and dd/dp may then be found efficiently by fitting a cubic spline to the 
values of L at each *&. 



For further details of the numerical integration of the holonomy equation ( |3.2| ) we refer 
the reader to [p3|l . In particular it is explained there how the currents Ri can be found 
directly from a holonomy equation involving the differentiated instanton data; the current 
R$ may be treated analogously. This is preferable to taking numerical derivatives of the 
Skyrme fields. Fields and currents are generated on a mesh of 50 x 50 x 50 points. However, 
because of the discrete symmetries Q4.2Q , it is sufficient to compute the fields and currents 
on a mesh of 25 x 25 x 25 points. 

Consider now the calculation of the potential, metric and static deuteron properties from 
the fields U(p, x) and the currents Ri(p, x) and R p (p, x). It is convenient to define quantities 
Mo, u, a,i and a p via 

U = Uq + IU-T 

Ri = icii-T 

R p = idp-T , (A. 6) 

and to introduce the following abbreviations: 

f(ct) = u x ol + UqOl 
F( a ,(3) = /(a)(u./3)-/(/3).(u.a). (A.7) 

Then the integrand of the potential energy (|2~T0|) takes the form 

cii-cLi + |((ara i )(a 3 --a :/ ) - (ai-aj)(ai-aj)) , (A. 8) 



and the terms that appear in the integrands of the metric coefficients ( |5.4|) - ( |5.7|) can be 
written as 

tr(x xR)l = -2& r &! 

tr{[T U (j\tf? = 8(ul + uJ) 
tr([r 3 ,U}U^xx R) 3 ) = 4if 3 (b 3 xu) 

tr([(x x R^ki^x x R)i,Ri]) = -8((bi x ai ) 2 + (xj + x 2 3 )(a 2 x a 3 ) 2 ) 
ti([[r u U}U\RMri,U]U\R t }) = 32 ^(a?(^ + u\) - ft{u x a t )) 

i 

tr([[T 3 ,U]fr,Ri][(xx R) 3 ,Ri]) = -mF 3 (a h b 3 x a<), (A.9) 



together with cyclic permutations. 
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The formula ( [2.12 ) for the baryon density B°(p,x) for fields in standard orientation is 
now 

B° = _(JL)a r 03X03, (A.10) 

and the static deuteron properties may be calculated from the following three moments of 
the baryon density: 

hip) = ljd 3 x x 2 B°(p, x), z = 1, 2, 3. (A.ll) 

Explicitly one finds 

r 2 = h + h + h 

rms - 1 - ' * ' u 

Qn = 2/ 1 -J 2 -J 3 

Q22 — 2J 2 — Ii — I3 

mil = iih+h) 

m 22 = |(/i + / 3 ). (A. 12) 

All these integrals have to be computed for several values of p. The integration over 
R 3 is facilitated by mapping the whole of R 3 bijectively to the finite cube C = [—1, 1] x 
[—1,1] x [—1,1]. The mesh used is uniform on C, but the mapping between R 3 and C 
changes with the Skyrmion separation. In keeping with the earlier discussion, consider 
coordinates (xi,X2,xs), with respect to which the Skyrmions are separated along the 1- 
axis. The cube C has coordinates (ci, c 2 , C3), each taking values in [—1, 1]. The form of the 
mapping is motivated by the observation that at large separations the Skyrmion positions 
are approximately (±Li,0,0) and their size is approximately y/L 2 L, where L\ and L 2 are 
the lengths of the semi-major and the semi-minor axes of the ellipse in the Harthorne data 
of figure 1. The aim is to concentrate mesh points around the Skyrmion centres, where the 
integrands make their largest contributions. Explicitly, 



(A.13) 

for j = 2, 3, where k is a dimensionless parameter that controls the degree to which points 
are concentrated near the Skyrmion centres; all calculations here had k = 0.4. The relation 
between c\ and x\ comes in two parts: 

x 1 = c 1 ^AL 1 {1-c 1 )-{k^L 2 ~L){1-2c 1 )^ for Qi 6 [0, |) (A.14) 



and 



L. + it^fa-i) for , 

[i-ciY 
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The two pieces fit together smoothly at c\ — ~. The expressions for c\ negative are chosen 
so that Xi(ci) is an odd function. 



After the mapping, the integrals are calculated via tricubic interpolation, as in [23]. The 



factors (1 — q) 2 in the denominators of (|A.13|) and ( A.15|) ensure that all the integrands have 



vanishing gradients at the boundary of C. (The denominator in [23|] was simply (1 — c, 
which is fine for integrands falling off at least as fast as r~ 6 for large r; here, some of the 
metric integrals fall off like r -4 , and this requires a different power if the tails are to be 
handled correctly.) 



We have computed the integrals ( |A.8|) , (|A.9|) , ([A. 10|) and ( |A.11|) for sixteen different values 



of p. The corresponding values for d are all the multiples of 2° in the range 2° < d < 30°, 
and an additional point at d = 1°, included to have a better check against the asymptotic 
expressions ( |5.19|) - ( |5.21|) . The range of values of p covered in this way is [1.71, 15.3] (the size 
of an isolated Skyrmion in these units is ~ 1). To compute all quantities for each separation 
takes approximately one hour on a workstation with an R4000 processor. There is good 
agreement asymptotically with ( |5.19| ) - ( |5.21| ) and the additional symmetries ( |5.15| ) and 



( [5.1 6| ) when p = po are also clearly present. The integrated baryon number density never 
leaves the range [2.0003,2.0049], and there is no evidence for regions of negative baryon 
density. 
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Figure Captions 



1. The Hartshorne data for instanton-generated two-Skyrmions in the attractive channel. 
The circle and the ellipse are shown in standard orientation; the isosceles triangle with 
vertices X 1 = (0,L,0,0),X 2 = (-L cos —L sin $,0,0) and X 3 = (L cos —L sin 0, 0) 
specifies the JNR data used in our calculations. 

2. Spatial moments of inertia. 

3. The spatial moment of inertia a 2 , shown on a smaller scale than in figure 2. 

4. Iso-spatial moments of inertia. 

5. The cross term w which couples spin and iso-spin. 

6. The radial metric coefficient f 2 . 

7. The potential V (solid line) and the effective potential V eB (dashed line) which occurs in 
the radial Schrodinger equation for the deuteron. Both are plotted in geometrical units. 

8. a) The radial part Ud of the deuteron wavefunction; for this plot the normalisation is 
chosen so that u(p ) = 1. 

8. b) The radial probability distribution u 2 abcABC f for the deuteron state, normalised so 
that dp ujabcABCf = 1. 

9. a) Equally spaced density contours for the spatial probability distribution p d of nuclear 
matter in the deuteron state ^d- The distribution is axially symmetric about the 3-axis, and 
shown here in the xiX3-plane. Both x\ and £3 are measured in fermi. 

9.b) Density contours for the spatial probability distribution of the nucleons calculated 
in a conventional potential model of the deuteron, with the nucleons treated as point-like 
particles. (From plf , with kind permission of Oxford University Press) 
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